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Abstract. Resummations of inhnite sets of higher-order perturbative contributions are often 
needed both in thermal held theory and at zero temperature. For instance, the behaviour of 
the Standard Model (SM) effective potential extrapolated to very high energies is known to 
be extremely sensitive to higher-order effects. The 2PI effective action provides a systematic 
approach to consistently perform such resummations. However, one of its major limitations 
was that its loopwise expansion introduces residual violations of possible global symmetries, 
thus giving rise to massive Goldstone bosons in the spontaneously broken phase of the theory. 
We review the recently developed symmetry-improved 2PI formalism for consistently encoding 
global symmetries in the 2PI approach, and discuss its satisfactory held-theoretical properties. 
We then apply the formalism to study the infrared divergences of the SM effective potential due 
to Goldstone bosons, which may affect the stability analyses of the SM. We present quantitative 
comparisons, for the scalar sector of the SM, with the approximate partial resummation 
procedure recently developed to address this problem, and show the quantitative discrepancy 
of the latter with the more complete 2PI approach, thus motivating further studies in this 
direction. 


1. Introduction 

In thermal field theory, finite-order perturbative expansions break down at high temperatures 
and one needs to resum higher-order contributions to deal with this problem. On the other hand, 
also at zero temperature there are situations where higher-order effects may potentially play an 
important role, even in a small-coupling regime. For instance, it has recently become well known 
that the behaviour of the Standard Model (SM) effective potential, extrapolated to very high 
energies, is extremely sensitive to the physics at the electroweak scale mw- Thus, a formalism 
to incorporate and resum higher-order effects in a consistent manner is highly desirable for of 
both thermal and non-thermal applications. 

A natural framework to address such problems is the formalism introduced by Cornwall, 
Jackiw and Tomboulis (CJT) |3j. Its simplest version is known as the Two-Particle-Irreducible 
(2PI) effective action. This is an effective action expressed not only in terms of the field, 
but also in terms of dressed propagators. When one considers a given truncation to the 
2PI effective action, at any given order of the loopwise expansion, the 2PI effective action 
resums, automatically, an infinite set of higher-order diagrams induced by partially resummed 
propagators, without the danger of over-counting graphs. 

So far, the 2PI formalism has been mainly used in thermal contexts, both at equilibrium 13 El 
HlEli and out of equilibrium uniiiiiiiiiis]- One of the major limitations for its application to 




zero-temperature problems like, for example, the study of the SM effective potential, has been 
the well-known difficulties of the CJT formalism to incorporate symmetries in a satisfactory 
way. In particular, global symmetries are not exactly maintained at a given loop order of the 
2PI expansion, since they get distorted by higher-order effects, and this results in violation of 
the Goldstone theorem [aiis] by higher-order terms, giving rise to a massive Goldstone boson 
in the Spontaneous Symmetry Breaking (SSB) phase of the theory [T6l [T71 HU [H] . In the past, 
several studies were presented in the literature, attempting to provide a satisfactory solution to 
this problem (see references in [20]). 

We have recently developed a new symmetry-improved CJT formalism |20| that addresses this 
long-standing problem and allows one to consistently use the 2PI effective action to study theories 
with global symmetries. In the symmetry-improved 2PI formalism, the Ward Identities (WIs) 
associated with global symmetries provide additional constraints for the extremal solutions of the 
fields and propagators. This formalism has a number of satisfactory field-theoretical properties, 
and actually provides an improvement over the standard GJT formulation, in the sense that the 
behaviour expected for the full theory is recovered already at a low orders of approximation. 
The symmetry-improved formalism has been also used to study the chiral phase transition |21j . 
and it has recently been extended to higher nPI effective actions [22| , confirming its consistency 
also when one goes beyond the 2PI approach. 

It has been recently pointed out [23| that the SM effective potential suffers from infrared (IR) 
divergences due to the Goldstone bosons of the electroweak gauge group, which appear at three- 
loop order. They appear at two loops for the derivative of the potential, i.e. for the minimization 
condition that fixes the value of the Higgs vacuum expectation value (VEV). This problem needs 
to be addressed for a two-fold reason. From the conceptual point of view, the effective potential 
should be well-defined for all values of the Higgs background field (j). On the other hand, 
from the quantitative point of view, the appearance of these IR divergences formally lowers the 
order of the involved three- and higher-loop contribntions, thus breaking down perturbation 
theory and potentially giving a significant contribution to the threshold corrections to the VEV. 
Because of the extreme sensitivity of the effective potential extrapolated at high energies to the 
matching conditions at the electroweak scale, this could result on quantitative sizable effects 
on the stability analyses of the SM potential. These issues have been recently addressed by 
means of an approximate partial resummation procedure, developed in [241125| . In addition to 
confirming the disappearance of the IR divergences within our approach, we will show that the 
more complete 2PI analysis, based on our symmetry-improved formalism, gives quantitatively 
different results, at least for the scalar sector of the SM. This suggests that even for the full 
SM there might be quantitative discrepancies between the approximate resummation scheme 
of [241 [25] and the 2PI approach, thus motivating further studies in this direction |26j . 

The layont of the paper is as follows. In Sections and [3 we review the standard and 
symmetry-improved 2PI formalisms, respectively. In Section ^ we consider the first non¬ 
trivial truncation of the effective action, known as the Hartree-Fock (HE) approximation, for 
a scalar 0(2) model, and demonstrate some of the satisfactory field-theoretical properties of 
the formalism. In Section we go beyond this approximation and show that the symmetry- 
improved approach describes correctly the thresholds of the Higgs and Goldstone self-energies, in 
agreement with nnitarity requirements. Moreover, we demonstrate explicitly how the formalism 
implicitly resums sets of processes at arbitrarily high order. Finally, in Section we use 
the symmetry-improved 2PI effective action to study the issue of IR divergences in a global 
SU(2) X U(l) scalar model. We draw our conclusions in Section]^ 


2. CJT Effective Action 

In this section we briefly review the GJT formalisms. For concreteness, we consider a 0{N) 
scalar model and show that, as opposed to the IPI formalism, loopwise truncations of the 
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CJT generating effective action lead to violation of the Goldstone theorem through higher-order 
terms. 

The model under consideration here is the O(A^) scalar model described by the Lagrangian 

m = i ^ (<pv - j (<pV , ( 2 . 1 ) 

where , (j)"^ , ■ ■ ■ , cj)^) represents the 0{N) scalar multiplet and summation over the 

repeated indices i, j = 1, 2,..., is implied. At zero temperature, T = 0, the model has a SSB 
phase for m? > 0, with breaking pattern: O(A^) —)• 0(A^ — 1). As a consequence of the Goldstone 
theorem [HIE], there are — 1 Goldstone bosons and one Higgs particle H. For the simple 
case N = 2, which we are considering here, the field components may be decomposed as 

= {^) + H , = <^2 ^ G , (2.2) 


where we have denoted field operators with a caret. Here, G is the Goldstone field, which is 
massless at the minimum of the potential, whereas H is the Higgs boson, which is in general 
massive. 

The CJT effective action is the generalization of the IPI one, where in addition to the local 
source Jx, multi-local sources are introduced. In its simplest version, the 2PI formalism, we 
consider one local and one bi-local source, i.e. Jx and Kxy Therefore, the connected generating 
functional W[J,K] is given by 


W[J,K] 


—i In / P(/)* exp 


+ Jlcpx + 



(2.3) 


where 5'[(/>] = £[</>] is the classical action. Here and in the following, repeated spacetime 

coordinates will denote integration with respect to these coordinates. The background fields <f)\. 
and the connected propagator A^j/ are obtained by single and double functional differentiation 
of W\J,K] with respect to the source J*: 


5W[J,K] 



. dW[J,K] 

dJp. 5JI 




xy 


(2.4) 


We also have 


6W[J,K] 

5K% 



(2.5) 


The 2PI effective action A] is obtained as the double Legendre transform of W[J,K] with 
respect to J and K: 


r[</.,A] = W[J,K] - Jlct^x- + cjp^x^y) . (2.6) 

This has the explicit form j^: 

T[(P,A] = S[cj)] - I Tr(lnA) + ^ Tv(a(°)-^a) - zT^^^) ^ ^2.7) 


where A® = (5^5[(/)]/(J(/)^ Jc))^) is the inverse tree-level propagator matrix and P^-^) is the 
diagrammatic series of all two- and higher-loop 2PI vacuum diagrams in which the dressed 
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Figure 1. Infinite set of perturbation-theory diagrams implicitly resummed, in the 2PI 
formalism, by the self-energy on the LHS. Thick lines represent the full propagator, thin lines 
the tree-level one. 


propagator matrix A is used, instead of the tree-level one. The equations of motions (EoMs) 
are obtained as 




- P — 


6A^y 2 ■ 


( 2 . 8 ) 


In the limit of vanishing external sources, the physical solution is obtained by extremizing the 
2PI effective action r[(/). A]. In particular, the EoM for the propagator has the form 


A-i = a(°)-![</)] + n[0,A] , 


(2.9) 


where !![(/>, A] is the IPI self-energy, in which the propagator lines are given by the dressed 
propagator matrix A. When a given truncation of r[i^, A] is explicitly considered, this EoM 
implicitly resums an infinite set of perturbation-theory Feynman graphs, as shown in Fig. 

By considering the 0{N) invariance of the 2PI effective action r[(/), A] we may find a WI 
valid in the 2PI formalism at the extremal point of the 2PI effective action [20]: 


6^mA] , ^ 6^mA] 


6GrSG. 


+ 2 


y 


SGySAgH 


A 


HH 


- A 


GG 


= 0 


( 2 . 10 ) 


When 0(A)-symmetric truncations of the 2PI effective action are considered, the approximate 
solutions and Ag satisfy (2.10), which is the WI of the 2PI formalism, but not, in general, 
the Goldstone theorem, which is a WI of the IPI formalism. Unlike the latter, the WI (2.10) 
does not protect the masslessness of the Goldstone boson. In fact, as it has explicitly been shown 
in mm I19j . the Goldstone boson G is not massless in the HE approximation, manifesting 
itself as a pole at 0 in the Goldstone-boson propagator A^^{k). We stress here that the 
exact solutions obtained from the complete 2PI effective action satisfy, instead, the Goldstone 
theorem, since the complete IPI and 2PI effective actions are fully equivalent at their extremal 
points Hj. 


3. Symmetry-improved 2PI Effective Action 

In this section we present the symmetry-improved formalism for the 2PI effective action P)^, A], 
which respects the Goldstone theorem for any 0(A)-symmetric truncation of A]. 

From the discussion in the previous section, it is clear that when a given truncation of the 2PI 
effective action rtr[</>, A] is considered, the solution obtained extremizing A] with respect 

to (j) and A fails to satisfy the Goldstone theorem 


A 


-1,GG 

xy 


= V A 


-1,GG 


(fe) 


k=0 


= 0 


(3.1) 


The symmetry-improved EoMs are obtained by imposing the IPI WI (3.1) directly as a constraint 


for the Goldstone propagator. This can be achieved, in a variational formulation, by performing 
a constrained extremization of the 2PI effective action, i.e. by finding the extremum of 


mA,£] = rt,[</.,A] - A-J’GG^ 


(3.2) 
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where is a Lagrange-multiplier field. In 
the 2PI EoM with the IPI one 


it is shown that this procedure results in replacing 


H. 


H 


= 0 


= 


5P^ 

(3.3) 


(3.4) 


The latter is explicitly given by the IPI WI 

dPipiid)] 

Therefore, assuming homogeneity, the symmetry-improved EoMs to solve in the SSB phase are 

drto[u,A] 


6A^i{k) 


k=0 


= 0 


= 0 , 


(3.5a) 

(3.5b) 


namely the standard 2PI EoM for the propagators and the Goldstone theorem itself. In the 


symmetric phase of the theory we have u = 0 and only (3.5a) needs to be solved. 


For homogeneous background field values (p away from the minimum v, i.e. for p ^ v, we 


may define a symmetry-improved effective potential Ves{(j)) by means of the IPI WI (3.4), seen 


as a differential equation defining the improved potential. We thus have the definition [20! 


dp 


(3.6) 


where A = 0;d>) is the Goldstone component of the solution of the EoM for the 

propagators (3.5a) for generic p ^ v. This definition is manifestly compatible with the EoMs at 


p = V (3.5). The integral form of (3.6) is 

Vesi^) = - [\pPA-^’^^{k = 0-p) + Eeff(d. = 0) 


(3.7) 


We show in |20j that the additive constant I4fr(d’ = 0) can be chosen such that the formalism 
is thermodynamically consistent in the sense discussed first by Baym in 1271. However, for the 
purposes of this work, we do not need to determine it. We refer the interested reader to [20] for 
a more detailed discussion. 


4. The Hartree-Fock Approximation 

In this section we apply the symmetry-improved approach outlined in the previous section 
to the HE approximation of the 2PI effective action. We show that the predicted Goldstone 
boson is massless, as it should be, and that the phase-transition is second order already in 
this approximation. These predictions are in agreement with general field-theoretic properties 
that hold for the full effective action of the theory, thus showing that the symmetry-improved 
formalism gives actually improved predictions as compared to the standard truncated effective 
action. 

The HE approximation consists in retaining only the double-bubble graphs (a)-(c) of Fig. 
The unrenormalized 2PI effective action in this approximation is given by 

r„p[.,A«,A=] = - Jf-) - (Tt(lnA'') - iTr(lnAG) 


+ 



—6iX 


— i 


iKx 


.H 




. —6fA 


i 


8 




G 

XX ’ 


(4.1) 
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Figure 2. Unrenormalized two-loop contributions to r[(;/!), A], with thick lines denoting dressed 
propagators. The HF approximation consists of the graphs (a), (b) and (c) and the sunset 
approximation includes the graphs (d) and (e). 


In the above, we have simplified the notation for the diagonal propagators as = A^^ and 

= A^g. 


4 . 1 . Renormalization 

It has been shown in |28 | l29 t l30] that the 2PI effective action is renormalizable with temperature- 
independent counterterms (CTs). 

To renormalize the theory, we consider all parameters occurring in the CJT effective 
action ( |4.1| ) to be bare (denoted with the subscript B): 

(I/b = , ml = Z-\m^ + 6m^) , As = Z-^{X + (5A) , A^ = ZA^^ . 

(4.2) 

In the 2PI formalism, there are two distinct 2-point operators appearing in the effective 
actions, corresponding to the two different derivatives 

S^rtr[cP,A] 6Ttr[<P,A] 

' 6A ’ ^ ’ 

When a generic truncation of r[(/). A] is considered, these two operators are independent 
(including their divergences), and need to be renormalized by two different CTs [3Uj . 
Analogously, in the truncated 2PI formalism we have the independent 4-point functions 


6^rtr[cl),A] 6^Ttr[<P,A] 

54)6(1)5(f)6(1) ’ 5(t)5(t)5A ’ (iA(5A 


Moreover, in the 0(A) model these operators appear in two different 0(A)-invariant 
combinations, denoted by A and B. It is shown in |30] that only the standard perturbation- 
theory coupling-constant CT 5X is required to renormalize the higher-loop graphs of the 2PI 
effective action, when going beyond the HF approximation. Hence, a finite number of CTs is 
needed to make the effective action finite, thus guaranteeing its renormalizability. With the 
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above considerations, the effective action (4.1) reads: 


rHF[D.A".A°| = 




-Tr 

2 


-Tr. 

2 


I Z9^ + ^3A + (5A^ + 2(5Af ^ A^| 

I Z+ (^\ +bXi'^v^ — + 5rn^ 


. -i (3A + <5A^ + 25A^) . 


— i 


zA^CiA^l. - ^ -2^(A + c^A^) 


"^XX '"^XX 


^xx '^‘^xx 


-i (3A + 6 X 2 + 2(5Af) Q Q 

i . LL^rf.rf. tL\rf 


•^xx '"^XX ’ 


(4.5) 


where we may set Z = 1 at this order of loop expansion. 

By canceling separately divergences and sub divergences proportional to temperature- 
dependent terms, in [20] we find the following CTs in the HF approximation, in the MS scheme 
of dimensional regularization d = 4 — 2e: 


3X^ = <5Af 


2A2 

IGvr^e 


6A 8A2 ’ 

IGir^e (lOvr^e)^ 


(4.6a) 


dAf = dAf 


2A2 1 


IGvr^e 2A ’ 
IGvr^e 


Smf 


4XrrP 1 
IGvr^e 4A 
IGtt^c 


(4.6b) 


The above T = 0 CTs are sufficient to renormalize the EoMs for the propagators, also when 
thermal effects are considered. 


4 . 2 . Thermal Phase Transition 

In the HF approximation, the self-energies are momentum independent. Therefore, we may 
parameterize the propagators as = (fc^ — -|- ie)~^, where the effective Higgs and 


Goldstone masses, Mfj and Mq, depend only on the temperature T. 

In the symmetric phase of the theory the constraint (3.5b) is automatically satisfied, and the 
renormalized EoMs (3.5a) take on the form 


Mjj 

Ml 


—m 


—m 


. Ml 
+ 3A ^ 


1 Ml 

-In —— 

167r2 2m2 


X Ml 
+ A ^ 


— InM 

IGtt^ 2m2 


X Ml 
+ X ^ 


X Ml 

-In —— 

IGtt^ 2m2 


X Ml 
+ 3A ^ 


1 Ml 

-In —^ 

IGtt^ 2m2 


-h 3A 


+ A 


n{uj 


H\ 


0 J^ 


H 


n{(jj- 


H\ 


UJ- 


H 


+ X 


n(a;j 


G\ 




G 


+ 3A 


n(a;^) 


w, 


G 


(4.7a) 

(4.7b) 


where = /d^k/(27r)^, cuk = Vk^ -|- is the on-shell energy of the (quasi)particle, and 
n(u;) = (e^/^ — 1)“^ is the Bose-Einstein distribution function. In the symmetric phase we 
obtain, as expected, a single solution with Ml = Mfj. 

In the HF approximation the constraint (3.5b) reads 


V mI = 0 . 


(4.8) 
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Figure 3. The values of Mfj, Mq and the VEV v, as functions of T, predicted in the HF 
approximation of the symmetry-improved 2PI formalism. 


In the SSB phase of the theory this implies Mq = 0, yielding the EoMs [20j 


Mfj = 3Au^ — + 2>\ 


Mjj , Ml 


IGvr^ 


In 


2r5 + 


n{u}^) 


UJ- 


H 


+ A 


n(w^) 


0 = Ai," - m2 + In^ + A / 

Idtt^ 2m^ uj- 


H 


+ 3A 


UJ- 

n{uj^) 


k 

G\ 


LJ, 


G 


(4.9a) 

(4.9b) 


We point out that have chosen the MS mass scale /i such that the tree-level relations MJj = 2m? 
v'^ = rn? jX are satisfied at T = 0. The mass-gap equations (4.9) are solved analytically to be 


, ,9 9 O AJ. 

Mjj = 2m - — , 

m2 = 0 , 

2 ^ ^ _ Mjj Mjj 

A 167r2 2m2 


[ n{uJ^) _ ^ 
k uj^ 12 


(4.10a) 

(4.10b) 

(4.10c) 


In Eig. we exhibit the dependence of the squared thermal masses, Mjj and Mq, and the 
thermally-corrected VEV v, as functions of the temperature T. The parameters of the model 
are chosen such that Mh = 125 GeV and v = 246 GeV at T = 0. We observe that we predict a 
second-order phase transition at T = Tc = y/3v{T = 0), already in the HF approximation, 
in agreement with theoretical expectations to all orders. This is in sharp contrast with 
the incorrect first-order phase-transition predicted in the HF approximation by the previous 
approaches [3I1E3E3]- 


4-3. Symmetry-improved Effective Potential 

Let us now calculate the effective potential at high temperatures in the HF approximation of the 
symmetry-improved formalism. Extending the renormalized EoMs (4.9) from u —>■ (/>, we obtain 
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Figure 4. Symmetry-improved HF effective potential in the high-temperature approximation. 


The large dots denote the minimum solutions cj) = v, obtained in Section 4.2 


MU^) = + 3A In 


y 2m^ J Ihvr^ y 2m^ J 


+ 3A 




'k 


+ A 


nluj 




MU4>) = + a 


+ A 


nluj 




(</>) 


2m? 


+ 3A 


f 

'k 


+ 3A^^ln 

167r2 


MU<P) 

2m^ 


(4.11a) 


(4.11b) 


1 

(j) d(j) 


= mI{4>) . 


(4.11c) 


The first two equations are the EoMs for the propagators, for general (j) ^ v. The latter 


equation (4.11c) results from the definition of the symmetry-improved effective potential (3.6), 


i.e. the IPI WI for (p ^ v. Figure presents our numerical estimates for the high-temperature 
symmetry-improved effective potential, as functions of p, for different temperatures T |20j . 
Again, a second-order phase transition is described, already in the HF approximation. The fact 
that Ve{i{(f)) acquires an imaginary part when p < v signals the instability of the homogeneous 
vacuum in this region 


5. Threshold Properties 

In this section, we study the threshold properties of the Higgs and Goldston particles, by 
including the contributions from the sunset diagrams (d) and (e) in Fig. In particular, 
we will show that the resummed Higgs- and Goldstone-boson propagators predicted within our 
symmetry-improved CJT formalism exhibit the correct threshold properties arising from on- 
shell Higgs and Goldstone particles in the loop. Therefore, we explicitly demonstrate that the 
symmetry-improved approach is consistent with the optical theorem and unitarity. 

One of the common approaches in the literature, when studying the issue of symmetries in 
the GJT formalism, is to define an additional 2-point function as 


A 


-1 

ext 


6(p 64> 


(5.1) 


9 
























H G H G 

^ Q + Q + -Q- + -Q- 

H G 

H G H 

G 


Figure 5. Higgs and Goldstone EoMs in the sunset approximation. 


where A{(f>) is the solution of the standard 2PI EoM for the propagator. This 2-point function is 
sometimes referred to as the external propagator. It can be shown that this function satisfies the 
Goldstone theorem, since the second derivative (in the Goldstone direction) of any symmetric 
(/>-only functional does so [20] . Thus, it is often claimed that this function should be considered 
as the true approximation to the propagator, in the GJT formalism, rather than A. However, 
the problem of this approach is that what appears in the diagrammatic series of r[</>. A] is the 
propagator A, not the external propagator Aext- The former does not satisfy the Goldstone 
theorem, in the standard GJT formalism, and therefore the thresholds of the particles are 
described incorrectly, since the Goldstone bosons propagating within quantum loops are massive. 
In this sense, this approach is not capable of describing the Goldstone bosons as consistently 
massless quantum-mechanically, i.e. within loops. From the above discussion, it is clear that 
studying the threshold properties described by the symmetry-improved 2PI effective action is an 
important step in establishing the consistency of the formalism and show further its advantages 
as compared to previous approaches. 

As diagrammatically represented in Fig. the symmetry-improved 2PI EoMs derived in the 
sunset approximation are given by 

A-^'^(p) = p^ - (3A + + 25Af)+ dmj - {3\ + 6\^ + 25 X^)Th 

- (A + aA) To + 1 + i , (5.2a) 

I 2 % 2 

A“^’^(p) = p^ — (A-I-dA^) + rn^ + 5m\ 

- (A + 5\^) Th - (3A + 6X^ + 25Af) To + - {-2iXvf 1hg{p) , (5.2b) 

i 

uA“^’'^(0) = 0. (5.2c) 

Here, we have abbreviated the loop integrals as follows: 

Ta = jf'' f iA“(/c) , Iab{p) = [ iA°‘{k - p)iA^{k) , (5.3) 

Jk Jk 


where a,b = H, G, Inp^ = In/r^ -|- 7 — ln( 47 r) and p is the MS renormalization scale. In (5.2) we 
have introduced the shorthand notation ^d^k/(2 'k)^. Details of the renormalization and 

the numerical approach to solve these self-consistent nonlinear equations are given in [20] . Here, 
we only outline the physical content of their solution. 

It is useful to define the effective energy-dependent squared masses M^^g(s) as 


A ^’^/^(s) = s - 


(5.4) 
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Figure 6. Numerical solutions for Re[M|f(s)] — (left frame) and \m[M‘fj{s)]/{T hMh) 
(right frame), where Mh and Th are the Higgs-boson pole mass and width, respectively. The 
dashed lines are the one-loop results in the IPI formalism. 



Figure 7. Numerical solutions for Re[Mg(s)] (left frame) and lm[MQ{s)]/{T hMh) (right 
frame). 


where s = is the Lorentz-invariant energy-squared parameter. In Fig. and Fig. we plot 
the dispersive (real) and absorptive (imaginary) Higgs- and Goldstone-boson mass squares, as 
functions of s. We see that there is a non-vanishing absorptive part ImM^(s) that results 
from the on-shell decay of the Higgs particle into two Goldstone bosons, i.e. H —)• GG. The 
threshold for this process is at s = 0, thus demonstrating that the Goldstone bosons in the 
loop are consistently treated as massless within our symmetry-improved 2PI formalism. This is, 
again, in sharp contrast with previous approaches in the literature, e.g. with the results found 
in m for the absorptive part of the external propagator, where the Goldstone boson exhibits a 
non-zero mass in the loop, as discussed above. 

The dashed lines in Figs. and show the predictions obtained at one-loop level in the IPI 
formalism. We point out that the kinematic opening of the thresholds is very sharp, as opposed 
to the smooth thresholds predicted by our symmetry-improved GJT formalism. As illustrated 
in Fig. this is due to the fact that the 2PI formalism automatically resums infinitely-many 
higher-order processes, without the need of explicitly considering them, including also processes 
that take place below the IPI threshold. 
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Figure 8. Some of the on-shell processes whose description is included in the absorptive part 
of the 1-loop CJT self-energy on the LHS. 


6. The Infrared Divergences of the Standard-Model Effective Potential 

In this section we describe the issue, recently pointed out in [23], of the IR divergences of the SM 
effective potential due to the Goldstone bosons of the electroweak group. Then, by considering 
the scalar sector of the SM, we show how the symmetry-improved 2PI effective action can be 
used to study this problem, and compare the results obtained in this framework with the existing 
approach OEni in the literature. More technical details of the 2PI analysis will be given in a 
forthcoming publication [26]. 

The effective potential 14ff of the SM, calculated in perturbation theory, suffers from IR 
divergences due to the appearance of Goldstone bosons in ring diagrams, as shown in Figure]^ 
At 3 loops the divergence is logarithmic, but it becomes more and more severe with increasing 
loop order. Moreover, these divergences are more severe when one considers the derivative of 
the effective potential For the latter, IR divergences start from 2-loop order (see Figure j^. 

In perturbation theory these divergences appear when the tree-level Goldstone propagators 
become massless, i.e. at the tree-level minimum of the potential. This is instead hnite, together 
with its derivatives, at the dressed minimum cj) = v. Nevertheless, this problem needs to be 
addressed, for the following reasons: 


(i) The effective potential 14fr(0) should be well-dehned for all values of </>, not only at its 
minimum (j) = v. Among other things, the off-shell 4> ^ v effective potential governs the 
dynamics of the background field in inflationary scenarios. 


(ii) At the dressed minimum cj) = v, the dressed masses of the Goldstone bosons vanish, and this 
implies that the tree-level mass thq is formally of the same order as the 1-loop Goldstone 
self-energy . Thus, since starting from 3-loop order the IR divergences of are of the 

form 1/(771^)”', n > 1 (see Figure 9) all these higher-loop contributions to are formally 
at 2-loop order. This means that perturbation theory breaks down and these diagrams can 
potentially have a significant impact on 2-loop results for i.e. on the state-of-the-art 
threshold corrections to the VEV v. In view of the extreme sensitivity of the SM effective 
potential, extrapolated at very high energies, to the matching conditions at the electroweak 
scale mmm. these issues can potentially affect the stability analyses of the SM. 


Moreover, a seemingly unrelated problem is that the tree-level mass of the Goldstone boson can 
be negative at the dressed minimum (p = v, thus generating an unphysical imaginary part for 
the SM effective potential at its minimum, which does not correspond to a true instability of the 
homogeneous vacuum. This fact should be contrasted with the discussion in Section [^ where 
we have shown that the symmetry-improved effective potential acquires an imaginary part only 
in the concave region corresponding to a physical instability. This suggests that a resummation 
of higher-loop diagrams is needed to address this conceptual issue as well. 

Finally, we point out that the location of the IR divergence depends on the value of the 
gauge-fixing parameter but the divergence is nonetheless present in any renormalizable 

gauge at the value of the field f for which rnQ{(p\ff} = 0. 
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Figure 9. IR behaviour of Goldstone-boson ring diagrams contributing to the effective potential 
and its derivative. The IR divergences start at three loops for already at two loops for 

its derivative dVes/dcf). 


6.1. Approximate Partial Resummation 

Let us outline the approximate partial resummation procedure presented in [2H [25] to address 
these IR issues. To facilitate an objective comparison with the 2PI approach discussed in this 
work, we will limit ourselves to the scalar sector of the SM, i.e. we will consider a global 
§U(2) X U(l) model with only the scalar Higgs doublet. Although this does not allow to draw 
quantitative conclusion for the complete SM case, this simplified model is sufficient to study 
this issue qualitatively and to compare the results obtained in the 2PI approach with the ones 
in 

The partial resummation procedure consists in considering ring diagrams, as in Figure]^ with 
insertions of 1-loop Goldstone self-energies He (A:). One approximates these 1-loop self-energies 
with their zero-momentum value nG'(O). With this important simplification, one can resum 
these diagrams by replacing the 1-loop Coleman-Weinberg contribution of the Goldstone bosons 
with 

3' 

2 

( 6 . 1 ) 

However, it can be shown that the derivative of this resummed term is still divergent. This 
problem can be solved by limiting nc(0) to contain only the terms not proportional to m^, i.e. 
by replacing 


K 


( 1 ) 


3 m; 


G 


eff,G 4(167r2) 




3_(m^+J^(0)) 

4 (IGvr^) 


2 r 


log 






UciO) Ug ^ nG(0) - ^^m^(log(m^/^2) _ _ (g_2) 

Notice that the subtracted term does not correspond to the contribution of a given diagram, 
but it is contained in the contribution of both the Goldstone tadpole integral and the Higgs- 
Goldstone sunset diagram. Finally, one needs to subtract from 14fr the diagrams that would be 
double-counted otherwise. In conclusion, one adds to the effective potential the term [2HI25] : 


■r r(resLim) 


3 (m^ -F Hg)^ 
4 (IGtt^) 


log 



3 

2 


K 


(d.c.) 


eS,G 


(6.3) 


where is the contribution of the double-counted diagrams, as discussed above. 

As we show in Figure |To| for the scalar sector of the SM, this procedure effectively resums the 
IR divergences present in dVef^/dej) at two- and three-loop orders. We include only the 3-loop 
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Figure 10. IR divergence of the derivative of the effective potential, as calculated in 
perturbation theory at 3-loop order, considering only the scalar sector of the SM. The gray 
dotted line is the tree-level contribution, the red dashed line the leading 3-loop one, whereas the 
black continuous line is the result of the approximate partial resummation procedure developed 

in [21 [25]. 


contributions coming from the Goldstone ring diagrams in Figure since they are proportional 
to l/m-Q and so responsible for the IR divergence. Moreover, as we discussed above, their order 
gets formally lowered at the dressed minimum. 


6.2. 2PI Approach to the Resummation of IR Divergences 

We now study these issues by means of the symmetry-improved 2PI formalism. This provides a 
more complete resummation, as compared to the approach outlined above, for several reasons. 
First, it is a first-principle approach, and no ad-hoc subtraction is needed. Second, as we 
are going to show below, the 2PI approach takes into account more topologies and does not 
necessitate to neglect the momentum dependence of the self-energy insertions that are resummed. 
Moreover, as we showed in the previous section, the threshold properties are correctly described 
within the symmetry-improved formalism. In particular, the Goldstone propagator in the 
formalism is massless at the dressed minimum of the potential and nevertheless the effective 
potential will be shown to be free of IR pathologies. 

We consider the global SU(2) x U(l) scalar model. By including, in the EoMs, the 1-loop (HF 
-|- sunset) GJT diagrams shown in the hrst line of Figure 11 we actually resum, automatically, 
a large class of diagrams, as depicted in Figure 12 Moreover, in order to be able to compare 


our results with the perturbative 2-loop calculation, we include the 2-loop GJT diagrams in the 
second line of Figure [m However, to simplify the treatment, we approximate the propagators 
appearing in these diagrams as the tree-level ones A((/)) ~ (</>). In this way, we take into 

account the full contribution of 2-loop topologies and, in addition, resum diagrams as the ones 
shown in Figure [Tsl Therefore, a much larger class of diagrams is included, as compared to the 
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Figure 11. Diagrammatic representation of the approximation scheme for the EoMs in 
Section [H 







Figure 12. Some of the topologies implicitly resummed by the 1-loop 2PI self-energies in the 
first line of Figure 



Figure 13. Some examples of the topologies resummed by including the 2-loop 2PI self-energies 
in the second line of Figure [TT| Notice that the propagators belonging to 2-loop 2PI topologies 
do not get dressed, because of the approximation A((^) ss used, for these, in Figure 
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method outlined in Section 6.1, and the momentum dependence of the resummed insertions is 
also retained. 

In the symmetry-improved 2PI approach, the IR divergences are absent by construction: 
IR divergences can be present only when two or more Goldstone propagators carry the same 
momentum, as in the ring diagrams shown in Figure In other words, the IR pathologies 
originate from chains of Goldstone lines with self-energy insertions between them. However, 
such topologies are necessarily 2-particle-reducible and thus do not appear in the diagrammatic 
series of r[(/i, A]. Therefore, the resummation of IR divergences is achieved automatically by the 
construction of the 2PI effective action. 

The EoMs can be easily obtained by generalizing (5.2) and including the 2-loop self-energies, 
in which A{(p) ~ A^*^^((/>). Their explicit form is given in Appendix A The numerical solution, 
in the vicinity of the dressed minimum (f> = v, is plotted in Figure 14 The black dots represent 


the numerical solution obtained in our approach. It is apparent that the results from the partial 
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Figure 14. The derivative of the effective potential, near its dressed minimum. We show results 
from 1-loop (green dotted line), 2-loop (blue dash-dotted line) and leading 3-loop (red dashed 
line) perturbation theory. The black continuous line is the result of the approximate partial 
resummation procedure discussed in Section 6.1 The black dots are the results obtained from 
the symmetry-improved 2PI effective action. 


resummation outlined in Section 6.1 do not reproduce the ones obtained in the more complete 
2PI resummation. We have checked explicitly that, expanding the EoMs at 2-loop order, we 
reproduce numerically the results coming from 2-loop perturbation theory. Notice also that the 
3-loop result is larger than what one would naively expect (roughly about A/167r^ times the 
2-loop one), because of the breakdown of perturbation theory, as discussed at the beginning of 
this section. In Figure 14 the symmetry-improved 2PI solution appears to be similar in size to 
the leading 3-loop result, potentially indicating that the effect of the resummation is small near 
the dressed minimum. However, this is a coincidence and depends on the class of topologies of 
graphs considered here. A more detailed study will be given in 


7. Conclusions 

The 2PI effective action provides a powerful theoretical tool to consistently resum inhnite series of 
perturbation-theory diagrams of different topologies. However, its loopwise expansion introduces 
residual violations of possible global symmetries by higher-order terms. In the case of global 
symmetries, this leads to the appearance of massive Goldstone bosons in the spontaneously 
broken phase of the theory. In this work we have reviewed the symmetry-improved CJT 
formalism, developed in [20] for consistently encoding global symmetries in loopwise expansions 
of the 2PI effective action. 

We have demonstrated, in a simple 0(2) scalar model, the key field-theoretical properties of 
the formalism. In detail, we have shown that the Goldstone bosons are described as massless, 
also within quantum loops, thus providing a consistent description of the threshold properties 
of the Higgs and Goldstone particles. Moreover, the thermal phase transition is at second order 


16 







already in the HF approximation. Thus, the behaviour expected for the full theory is recovered 
already in the first non-trivial approximation, contrary to other approaches in the literature. 

Thanks to the satisfactory field-theoretical properties of the formalism, we have applied the 
symmetry-improved 2PI effective action to the study of the IR divergences of the SM effective 
potential due to the electroweak Goldstone bosons. By limiting ourselves to the scalar sector 
of the electroweak group, we have confirmed that these IR divergences are indeed an artifact 
of perturbation theory, and are absent in the symmetry-improved 2PI formalism, as it should 
be. We have compared quantitatively our results with the predictions given by the existing 
approach in the literature, consisting in an approximate partial resummation of Goldstone- 
boson ring diagrams. We have shown that the latter, whilst correctly resumming qualitatively 
the IR divergences, does not reproduce quantitatively, at least in the simplified scalar model 
considered here, the results of the more complete 2PI resummation, near the dressed minimum 
of the effective potential. 

In view of the extreme sensitivity of the SM effective potential, extrapolated to very high 
energies, to the matching conditions at the electroweak scale, it is important to assess these 
issues in a more realistic way, going beyond the scalar sector of the SM. A detailed study of this 
matter, including the quantitatively most important contributions, will be given in a forthcoming 
publication. 
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Appendix A. Equations of Motion for the SU(2) x U(l) Scalar Model 

In this appendix we give the renormalized EoMs used in Section |6.2| to study the IR divergences 
of the effective potential. Since the scalar SU(2)i x U(l) model automatically possesses the 
custodial symmetry SU(2)i x SU(2)ij ~ 0(4), the HF and sunset contributions to the EoMs are 
easily inferred from the 0(2) case (see [20]). We find, in Euclidean space: 


A-i’^(p;(/.) 
A“i’^(p; 0) 

IdVeM 
(p d(j) 


= + 3A</,2 - + 3A7F" + 3A73“ - 18- QX^ 

+ , (A.la) 

= + X(p‘^ — m? -\- XTJl^ + hXT^^ — ^X^(p"^ZQfj{p) + Hg > 

(A.lb) 

= A“^’^(0;(/>), (A.lc) 


where the MS renormalized sunset and tadpole integrals are given by [20) 
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IGtt^ 
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M2 + 


Vg X^(p 
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B{p] 


(A.3) 


17 











where vh = 24, vq = 4 and we have introduced 


= m^{p) , with A^{k) = 


’ 


, = 1 - it ( -1.34391 

' \/3 I 2 36 i 


(A.4) 

(A.5) 


Here, we seize the opportunity to eliminate an error in the value of r] reported in 


The terms 


in the second line of (A.3) are the MS finite part of the CT integrals used to renormalize Ta by 

Their subtraction is needed to guarantee that the renormalization 


the method discussed in 

scheme used here matches the standard MS one in perturbation theory. Finally, the 2-loop 2PI 

2PI (2) 

self-energies H^:^ ^ (p; (f>) are calculated in perturbation theory by standard techniques 


2PI ( 2 '] 

We approximate them by their zero-momentum value {4>)^ since the error introduced 

in this way is expected to be negligible for the purposes of this work. Adopting the compact 
notation used in [35], in terms of the functions defined there we find 


- 6X‘^I{H,H,H) - 6X‘^I{H,G,G) 

- 216X^(j)‘^ - 72X^I{H',G,G) - 24AV ^(G', G,//) 

- 648AV /(//', F', FT) - 144AV^(^',G',G) - 24AV/(G', G', iF) , 

(A.6a) 

(167r2)2ng = 8X^(f>‘^ B{G,Hf - 24X^I{H,H,H) + 22X^I{G,H,H) 

- 16A2/(G,G,iF) + 6A2/(G,G,G) , (A.6b) 


where all functions have to be evaluated at zero momentum. 
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